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ABSTRACT 

We show that by analyzing the spectra of lensed broad absorption line quasars (BALQSOs), it is possible to 
reveal key properties of the outflowing gas in the inner regions of these objects. This results from the fact 
that each image of the quasar corresponds to a different line of sight through the outflow. This combined with 
dynamical estimates for the location of the flow, adds new information concerning the lateral, non line of sight 
structure of the absorbing gas. Here we consider a sample of lensed BALQSOs and note that the similarity 
of BAL profiles of different images of the same quasar implies that the flow is relatively isotropic on small 
scales. We show that its geometry is inconsistent with the ballistically accelerated spherical cloud model, and 
that wind models provide a better description of the flow structure. Furthermore, observations seem to disagree 
with naive interpretations of recent time-dependent wind simulations. This may hint on several important 
physical processes that govern the structure and dynamics of such flows. Future prospects for the study of 
quasar outflows with the effect of strong gravitational lensing are discussed. 

Subject headings: gravitational lensing — ISM: jets and outflows — galaxies: active — galaxies: nuclei — 
quasars: absorption lines 



1. INTRODUCTION 

The standard model for quasars consists of a central contin- 
uum emitting region (possibly an accretion disk), broad and 
narrow line regions (BLR and NLR respectively), and a torus. 
Among these components only the NLR can be resolved in 
nearby quasars to reveal a stratified structure and a cone-shape 
geometry. The physical nature of the other components is sub- 
ject to an ongoing debate, e.g., it is not known whether the 
BLR is in the form of virialized clouds (e.g., Kaspi & Netzer 
1999) or a continuously outflowing wind (e.g., Murray, Chi- 
ang, Grossman, & Voit 1995; hereafter MCGV). 

Observations indicate that outflowing gas is a common phe- 
nomenon in quasars. Such gas is intrinsic to the quasar, as 
has been shown for the high velocity (> 10^^ km s^'), broad 
absorption fine (BAL) flows in BALQSOs (Weymann 1977). 
The quasar unification scheme (e.g., Elvis 2000) suggests that 
BAL flows are ubiquitous and their detection (in some 15% 
of type-I objects; Hewett & Foltz 2003) requires a fortuitous 
line of sight to the central engine. 

There have been several unsuccessful attempts to constrain 
the location of BAL flows via the detection of variations in 
velocity or ionization, de Kool et al. (2001) have detected 
absorption from excited levels of iron which allowed location 
estimates in a subclass of BALQSOs (FeLoBALQSOs). An 
indirect measure of the gas location relies on the observation- 
ally motivated assumption that it is accelerated by radiation 
pressure force (Arav 1996). This has been the focus of many 
works concerning the dynamical modeling of such flows. 

Dynamical models consist of cloud (e.g., Chelouche & Net- 
zer 2001, hereafter CNOl), wind (MCGV, Proga, Stone, & 
Kallman 2000; hereafter PSK), and cloud-wind models (e.g., 
Arav, Li, & Begelman 1994; hereafter ALB). Some versions 
of such models can be brought to some agreement with obser- 
vations. They do not however agree among themselves con- 
cerning the mass flow rate and other important parameters. 
Thus, despite some 30 years of research, several key proper- 
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ties of BAL flows remain unknown. 

In this Letter we show that strong lensing can add invalu- 
able information on the structure of BAL flows: In §2 we de- 
fine our sample of lensed BALQSOs. §3 discusses the obser- 
vational constraints on several flow models. Future prospects 
and caveats are discussed in §4. Summary follows in §5. 

2. LENSED BALQSOS: GEOMETRY AND SAMPLE 

In gravitationally lensed systems, each image corresponds 
to a different line of sight to the central object. For typical lens 
configurations (figure^, the angle between the line of sights, 
A^, is of order of the separation angle between the images 
(A0; see e.g., Schneider, Ehlers, & Falco 2000). 

We consider 7 lensed and one lens candidate (UM425; see 
however Aldcroft & Green 2003 who favour a binary quasar 
interpretation for this system) BALQSOs whose properties 
were taken from CASTLES^ (Munoz et al. 1998) and Michal- 
itsianos et al. (1997), and are given in table 1. Assuming a 
typical type-I quasar continuum with ttox = 1 .4 (Green et al. 
1995) and a X-ray (0.1-lOkeV) continuum slope, = 0.9, 
the mean demagnified bolometric luminosity of the quasars 
in our sample is ^ 10'*^ erg s^' (Chartas 2000; this value de- 
pends on the assumed lens model and is somewhat insecure 
for double imaged lenses). 

Spectroscopic observations show that in most lensed 
BALQSOs, the absorption troughs are similar between differ- 
ent images to within S/N and resolution limits. This is prob- 
ably not a selection effect (due to the search criteria used for 
finding lensed objects) since the fraction of BALQSOs among 
lensed quasars is larger (Chartas 2000; or similar but not 
smaller) than their fraction among quasars (Hewett & Foltz 
2003). This means that BAL flows are relatively isotropic 
on small angular scales. Nevertheless, for two objects in our 
sample (APM08279H-5255 and the lens candidate UM425) 
there seem to be noticeable differences in the BAL troughs 
between the images. We discuss the implications of these re- 
sults on various dynamical models in §3. 

^ http://cfa-www.harvard.edu/castles/ 
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Fig. 1. — A schematic (not to scale) view of a double image lens where 
each line of sight to the continuum source traverses different sections of the 
BAL flow (see text). 

3. IMPLICATIONS FOR QUASAR OUTFLOW MODELS 

Here we show that by comparing BAL features of different 
images one can put constraints on existing flow models. This 
resuhs from the different characteristic angular scale of the 
flow predicted by different models. Below we discuss each 
class of models separately. 

3.1. Single cloud models 

This model was initially invoked to explain the BLR gas 
and only later applied to BAL flows. Several works (e.g., 
Blumenthal & Mathews 1979, hereafter BM79) have demon- 
strated a qualitative agreement between model and observa- 
tions in terms of the flow velocity and spectral features. 

In these models the clouds are spherical to a good approx- 
imation (BM79) and are confined by some external pressure, 
P (P oc is typically assumed where r is the distance from 
the central object; ALB). The clouds are photoionized by the 
quasar radiation field and their ionization structure depends 
primarily on the ionization parameter, Ux (defined as the ra- 
tio of photon density in the range 0.1-lOkeV to gas number 
density; CNOl). The terminal velocity of the cloud (vioooo; in 
units of lO'* km s^') is given by 

vioooo-0.1L4fri7'/'Mi/2, (i) 

where = r/lO'^ cm is the launching distance, L45 the bolo- 
metric luminosity in units of 10"*^ erg s^', and M the "force 
multiplier" (the ratio of the total to the Compton radiation 
pressure force; M{r) = const, for the assumed pressure pro- 
file; CNOl). For 10"^ <Ux < 0.1 and the typical velocities 
and column densities of BAL flows, M can be approximated 
up to a factor ^ 3 by 

M~0.3Zt/_,r', (2) 

where Z is the gas metalicity normalized to solar values. 

High S/N spectroscopy combined with detailed photoion- 
ization calculations can constrain the ionization parameter, Ux 
and the column density Nii (in units of 10^' cm^^; assuming 
Z = 1) of the outflowing gas. These seem to cluster around 
t/, - 10""^, A^2i > 1 for BAL flows (Arav et al. 2001) and 
Ux - 10-^ N21 > 1 for FeLoBAL flows (de Kool et al. 2002). 
Therefore, high ionization BAL flows that are launched from 
the centers (rn > 50) of luminous quasars (L45 = 10) would 
reach a terminal velocity < 15,000 km s^' (equationsn&|2j 
which is in qualitative agreement with observations. 

Combining equations[2&|2lwe obtain, 

ri7 ~ 3 X IQ-\^^qoqL45ZUx-\ (3) 

where vioooo is inferred from the blue wing of BAL troughs. 
This estimate for rn can be compared to those derived for 
FeLoBALQSOs: For FBQS 1044H-3656 (L45 = 10 vioooo - 
0.5) de Kool et al. (2002) conclude that rn ~ 2 x 10"* (see 
however Everett, Konigl, & Arav 2002) which is consistent 
with the dynamical estimate for this object (rn ^ 10^). 



TABLE 1 

Properties of the Lensed BALQSOs Sample 
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APM 08279+5255 


3.91 
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a,b 


SBS 1520+530 


1.86 


0.72 


1.56" 
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2 
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HE 2 149-2745 


2.03 


0.49 


1.7" 


103 


0.8 


c,e 


RX J 09 11. 4+0551 


2.8 


0.77 


0.48" 
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0.3 


c,e 


PG 1115+080 


1.72 


0.37 


2.3" 


25 


0.6 


c,f 


UM 425 


1.46 


0.6 


6.5" 




1.2 


g.h.i 


H 1413+117 


2.55 




1.35" 




0.9 


cj 


J 1004+1229 


2.66 


0.95 


1.54" 




1 


k 



The sample of lensed (and the lens candidate UM 425) BALQSOs: [1] Object's 
name, [2] Source redshift, [3] Lens redshift, [4] Image separation (arcsec- 
onds), [5] Time lag between images (days), [6] Maximum flow velocity 
(10000 kms^') inferred from the absorption line troughs, [7] References: 
(a)Chartas et al. 2002, (b)twin et al. (1998), (c)Chartas (2000), (d)Burud' 
et al. (2002), (e)Bade et al. (1997), (f) Michalitsianos et al. 1996, 
(g) Michalitsianos et al. 1997, (h) Green et al. (2001), (i) Aldcroft & Green 
(2003) ,(j) Hutsemekers (1993), (k)Lacy et al. (2002) 

It follows from the definition of Ux that the angle subtended 
by a spherical cloud with a column density Z^ ' A^2i (e-g-. Arav 
et al. 2001) is 

A(pc ~ 10UxrnL^^Z^^N2i arcseconds. (4) 
By combining equations|3]and|3we obtain 

A(pc ~ 3 X 10"^Vj,jqqqA^21 arcseconds. (5) 

Thus, for all objects in our sample the cloud's lateral size, 
~ l0^^rnA(pc cm, is smaller than the Schwartzschild radius 
(Rs) of the central black hole (BH; using the luminosity-BH 
mass relation of Kaspi et al. 2000). The fact that BAL flows 
only partially cover the ionizing source (whose size is ^ IQRs 
for a standard accretion disk; Arav et al. 2001) corroborates 
our size estimates. 

Inspection of table 1 reveals that for all objects, A^^ ^ AO. 
For clouds to show identical absorption features in all quasar 
images, their aspect (lateral to radial) ratio should reach 10^ 
for the more extreme cases in our sample (the cloud model 
is however consistent with the different BAL troughs in e.g., 
UM425; Aldcroft & Green 2003). Such oblate clouds would 
probably be destroyed over dynamical timescales and a coher- 
ent acceleration is therefore unlikely (BM79). 

Thus far, we have considered single clouds. Neverthe- 
less, the flow column density may be distributed among many 
clouds (this is likely given the smooth nature of BAL troughs). 
For simplicity we assume that the flow consists of n identical 
clouds, each with a column Z 'A^21 /«■ In this case, A^^. is re- 
duced by a factor n. This model is sometimes referred to as the 
tube model where the aspect ratio of the flow is much smaller 
than unity (contrary to what is implied by lensed quasar obser- 
vations). The probability for a perfect alignment of the tube 
to several lines of sight is negligibly small . Thus, more likely 
geometries are those of a sheet or a cone which are naturally 
accounted for by wind models. 

3.2. Clumpy wind models 

In this model (ALB), the flow is assumed to consist of nu- 
merous small clouds ("cloudlets"). Isotropy is assumed in the 
statistical sense, i.e., each line of sight corresponds to a dif- 
ferent realization of the same cloud distribution. Thus, each 
image of the quasar may show different absorption features 
due to statistical fluctuations in the number of clouds between 
different lines of sight which follow the Poissonian distribu- 
tion. Following ALB, we assume identical cloudlets. Denot- 
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ing the total optical depth by T and the fluctuations in the opti- 
cal depth by At, one can estimate (assuming 1 a fluctuations) 
the average number of clouds along the lines of sight, 

n ~ t^At"^ (6) 

for At < T. As the spectral signature of a single cloud can- 
not extend considerably beyond one thermal width in velocity 
space (typically ~ 10 km s^') it allows us to constrain, or put 
lower limits to the number of clouds in each velocity bin for 
non-saturated troughs. In most cases, BAL profiles of differ- 
ent images are identical within S/N hmits implying that the 
number of clouds in each velocity bin is » 1 and their total 
number 100. This information is orthogonal to that con- 
cerning the filling factor of flows. Furthermore, observing 
systematic differences over a wide velocity range would indi- 
cate changes in the global flow properties (in the cloud distri- 
bution function). Such variations seem to be impUed by the 
dissimilarities in BAL troughs for two objects in our sample 
(APM08279H-5255 and UM425; see §4.1) although higher 
quality data are needed to confirm this point. 

3.3. Continuous wind models 

Time independent continuous wind models usually assume 
gas which expands according to the spherical continuity con- 
dition p <x. Q,r~^v~^ (p, V being the gas density and ve- 
locity respectively, and fi the opening angle; e.g., MCGV). 
This assumption is justified for Q. < 4-7t. However, flows are 
likely to evaporate quickly through their rims for consider- 
ably smaller Q,, since the expansion timescale (o= y/^m/vs, 
Vj ~ 20 km s~^ is the sound speed) is smaUer than the dynam- 
ical timescale (oc rn/vioooo)- Put differently, for ^/Q,/4k 4C 
400v^QQQQ arcseconds, the gas would evaporate before reach- 
ing its terminal velocity. Therefore, a confining pressure must 
be invoked, and it is highly unlikely that the resulting continu- 
ity condition would mimic a spherical expansion. If the differ- 
ences in BAL troughs in APM 08279-H5255 and UM 425 are 
real, these imply that ^Q./4n ~ AO (otherwise the probabil- 
ity for detecting such differences would be negligible). Thus, 
at least for these objects (see also Chelouche & Netzer 2003), 
the assumption of spherically expanding flows is unwarranted 
and alternative geometries with very different spectral predic- 
tions are more appropriate [e.g., confined (Konigl & Kartje 
1994) and freely expanding flows (Chelouche 2003)]. 

3.4. Time dependent wind models 

Time-dependent calculations (PSK) show that Rayleigh- 
Taylor, Kelvin-Helmholtz, and radiative instabilities result in 
a filamentary sub- structure of the flow. Currently, there are 
no spectral predictions for their model and direct comparison 
with observations is complicated. Nevertheless, the results 
for the angular distribution of Ux and density (see figure 4a,b 
in PSK for the angular range subtended by high velocity gas) 
imply that the extreme (> 10 orders of magnitude!) variations 
over < 10° result from opacity rather than density effects. 
Thus, Ux is a measure of the transmitted spectrum. Assuming 
that Ux is a continuous function of angle down to small scales, 
then interpolation of PSK results predicts it should differ by 
many orders of magnitude over ~ 1 arcsecond scales. Such ex- 
treme variations would result in different spectral signatures 
between quasar images which are not observed. A more de- 
tailed study of this effect (e.g., by accounting for the finite size 
of the continuum source and the stratified flow structure) is 



limited by the low resolution of current simulations. This ap- 
parent contradiction between model and observations may be 
resolved if there is a typical scale (larger than rA8 ~ lO^'* cm) 
below which variations are suppressed (e.g., due to the dis- 
sipation of turbulence below some scale or by the stablizing 
effect of magnetic fields; e.g., Proga 2003). 

4. FUTURE PROSPECTS 

We have demonstrated a new method for investigating the 
structure of BAL flows. Similar reasoning can be applied 
to additional types of flows namely, narrow absorption line 
(NAL) and X-ray BAL (XBAL) flows. For NAL flows N21 ^ 
2 (e.g., Kraemer et al. 2002; but see also Hamaim, Netzer, & 
Shields 2000 for lower values) and vioooo ^01 (Crenshaw et 
al. 1999) and, therefore, A^^ > 2 arcseconds, i.e., the single 
cloud model cannot be firmly ruled out for typical lens sepa- 
rations. This stems from the fact that such gas lies at larger 
distances and is more dilute. XBALs have so far been ob- 
served in two objects (APM08279-H5155 and PG 1115-h080; 
Chartas et al. 2002, 2003) and seem to have N21 = 10^, and 
^10000 = 10. This gives A(pc ^ 10^^ arcsecond. Nevertheless, 
due to the complications associated with X-ray spectroscopy, 
it is difficult to compare XBALs of different images. 

Given several lines of sight, it may be possible to detect 
whether the flow is collimated and at what direction. From 
pure geometrical considerations, a collimated flow will show 
a velocity difference 5v ~ v5(p between two lines of sight 
5(p apart. For SBS 1520H-530, 5v ~ 0.1 kms-\ i.e., within 
the resolution of next generation spectrometers. This, com- 
bined with an independent measure for the orientation of the 
quasar (e.g., via the detection of a jet), may allow us to as- 
certain whether BAL flows are predominantly polar or planar 
and how they relate to other components of the quasar. 

Outflows are Ukely to rotate around the central object (e.g., 
MCGV) and it may be possible to detect spectral absorp- 
tion features sweeping across our line of sight over Keplerian 
timescales, 

TKep. ^ 10-Sroooo^45Z^/^f/7^/^ri/2A(pdays, (7) 

(F is the luminosity in units of the Eddington luminosity). 
Typically TKep. ~ 100 days for BAL flows in quasars with 
F = 0.1. Detecting such events would set more stringent lim- 
its on the geometry of the outflowing gas and may reveal the 
angular momentum balance in the inner quasar regions (per- 
haps even in the accretion disk). Such studies may also pro- 
vide better estimates for Hq since absorption events do not 
depend on the global geometry of the system and are less 
smeared in time compared to emission events. 

4.1. Caveats 

We have assumed that flow properties do not change on 
timescales shorter than the time lag, fiag, between quasar 
images. This requires that both the dynamical timescale, 
fdyn and the intrinsic flux variabihty timescale, fflux, are 
longer than ?iag. Indeed, the dynamical timescale, t^yn ~ 

v^^qqqL4sZUx^^ days is much longer than fj^g (table 1) for UV 
flows (fdyn <C flag for XBALs) and recent variability studies 
of high redshift quasars indicate that fflux ^ flag (Kaspi et al. 
2003). Thus, any differences between UV BAL troughs are 
probably not related to time-lag effects. 

An additional complication arises from micro-lensing (ML) 
by stars in the lensing galaxy which can differentially (de- 
)magnify the continuum emitting region (e.g., Wambsganss 
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1992). Its effect, being line of sight dependent (e.g., Ofek & 
Maoz 2003), may alter the shape of non-black troughs and 
result in apparent discrepancies between them. A more quan- 
titative discussion requires further study (Chelouche 2003, in 
preparation) and is beyond the scope of this Letter. 

5. CONCLUSIONS AND SUMMARY 

In this Letter we show that by analyzing the spectral ab- 
sorption features of spatially resolved lensed BALQSOs, one 
can obtain invaluable information concerning the geometry 
and sub-structure of BAL flows. Given the currently avail- 
able spectroscopic data we conclude that such flows are rela- 
tively isotropic on a few arcsecond scales. Our analysis sug- 
gests that the coherently accelerated cloud model (and the 
related tube model) are inconsistent with observations, and 
that BAL flows are probably made of one or several outflow- 
ing sheets or cones with lateral dimensions much larger than 
the typical separation between quasar images. The similarity 



of the absorption features between different images is incon- 
sistent with naive interpretations of time-dependent dynami- 
cal simulations. This may indicate a typical length scale be- 
low which sub-structure is suppressed. Future high resolution 
spectroscopy of lensed BALQSOs may reveal whether such 
flows are collimated and distinguish between polar and planar 
flow configurations. Multi-epoch observations may put fur- 
ther constraints on the lateral geometry and dynamics of BAL 
flows as well as lead to better estimates for Hq. Such infor- 
mation cannot be obtained from observations of non-lensed 
BALQSOs, and is imperative to our understanding of quasar 
structure and the physical processes near black-holes. 
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